Cluster coherent potential approximation for electronic structure of disordered alloys 
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We extend the single-site coherent potential approximation (CPA) to include the effects of non-local 
disorder correlations (alloy short-range order) on the electronic structure of random alloy systems. 
This is achieved by mapping the original Anderson disorder problem to that of a selfconsistently 
embedded cluster. This cluster problem is then solved using the equations of motion technique. 
The CPA is recovered for cluster size Nc = 1, and the disorder averaged density-of-states (DOS) 
is always positive definite. Various new features, compared to those observed in CPA, and related 
to repeated scattering on pairs of sites, reflecting the effect of SRO are clearly visible in the DOS. 
It is explicitly shown that the cluster-CPA method always yields positive-definite DOS. Anderson 
localization effects have been investigated within this approach. In general, we find that Anderson 
localization sets in before band splitting occurs, and that increasing partial order drives a continuous 
transition from an Anderson insulator to an incoherent metal. 



PACS numbers: 71.30.+h 7I.10.Fd 72.15.Rn 
I. INTRODUCTION 

The problem of atomic short-range order (SRO) and 
its effect on the character of electronic dynamics has been 
studied for many years now. It is relevant to the detailed 
understanding of the conditions under which a transition 
from metallic to an Anderson localized (AL), disordered 
insulator occurs with increasing disorder strength [1]. In 
the d = oo limit, the well-known coherent potential ap- 
proximation (CPA) [2] provides the exact solution of this 
Anderson disorder problem. However, by construction, 
CPA cannot access specific effects (quantum interference 
from short-ranged scattering potentials) leading to AL 
behavior. Such an endeavor demands explicit incorpora- 
tion of the dynamical eSect of SRO on carrier propagation 
in a disordered system, and thus, a cluster generalization 
of CPA. Such attempts have indeed been carried out [3,4] , 
but are extremely cumbersome numerically. Addition- 
ally, they do not always guarantee the correct analytic- 
ity properties of the Green's functions: the few which do 
succeed in this respect, like the travelling cluster approxi- 
mation (TCA) [4] and the dynamical cluster approxima- 
tion (DCA) [5] are extremely tedious technically. Fur- 
ther, the study of the effects of short-range order (SRO) 
on carrier dynamics has, to our best knowledge, never 
been attempted using these approaches. It is obvious 
that attempting to tackle the harder problem of atomic 
(or doping induced) SRO along with arbitrarily strong 
local electronic correlations, not to mention important 
aspects like multi-orbital character of realistic transition 
metal-based oxide systems, with these approaches would 
be extremely numerically time-consuming. Given this, it 
is imperative to develop semi-analytical routes as far as 
possible, resulting in much better (and easier) numerical 
tractability. 

On the other hand, many important results have been 
gleaned from field-theoretic studies of the Anderson tran- 



sition. Perturbative renormalization group (RG) [6] 
approaches and extensions thereof, work in the weak 
disorder regime, and are a priori inapplicable in the 
non-perturbative regime where the Anderson-Mott MIT 
would be expected to occur in the 3d case. As a re- 
sult, well-defined precursors of the MIT are observed even 
at very high temperatures, as experimentally demon- 
strated [7] in many systems in the form of scaling behav- 
ior of various quantites, breakdown of Matthiessen rule 
and the Mooij correlation. Thus, while the perturbative 
{d + e) [6] approaches have indeed provided wealth of in- 
formation, such approaches are insufficient at strong cou- 
pling, which is precisely the regime of interest for doped 
TMO systems, as well as systems like strongly interacting 
two-dimensional electron systems (in Si MOSFETs [8]) 
which have been found to undergo insulator-metal tran- 
sitions. This is because one is always effectively in the 
strong disorder regime in strongly correlated electronic 
systems, where the renormalized one-elecron band width 
is very small (caused by Hubbard band-narrowing) in the 
correlated metal (or Mott insulating states, where the 
band splitting a la Mott-Hubbard mechanism occurs). 
Consideration of such cases is out of bounds with per- 
turbative approaches, and this requires development of 
genuinely non-perturbative approaches which should be 
capable of: 

(i) extending the CPA to access Anderson localization 
effects, and, 

(m) having sufficient flexibility to incorporate effects of 
Mott-Hubbard physics via dynamical mean-field (DMF) 
or cluster-DMF approaches. 

In this paper, we devise a new cluster-CPA tech- 
nique that satisfies the above requirements. It is ex- 
tremely simple to implement (requires only the solu- 
tion of Nc coupled, non-linear equations for the Green's 
function for a cluster with Nc sites), captures the in- 
tracluster correlations exactly, and is suited to further 
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improvements (larger cluster size, incorporation of Mott- 
Hubbard physics). As we will show, it also reproduces the 
exact CPA limit [2] for the single-site cluster {d = oo). 

The CPA is tailor-made to describe electron dynam- 
ics in a system with purely diagonal disorder when cor- 
relations between spatially separated disorder scatterers 
can be ignored, a situation which is formally exact in 
d = 00. Extensions to include off-diagonal (hopping) dis- 
order within the CPA framework have been proposed [3] 
by a variety of authors. Here, we will propose a different 
(related to non-local CPA) Green's function technique. 
Using the EOM technique, we explicitly include dynam- 
ical effects of arbitrarily strong scattering from short- 
range correlated disorder potentials. In the process, we 
will make explicit contact with the problem of Anderson 
localization in disordered systems. 

II. MODEL AND SOLUTION 

The first step is the construction of a suitable cluster 
model Hamiltonian incorporating diagonal disorder. Mo- 
tivated by results of Ref. [9] , we generalize the Anderson 
disorder problem to finite dimensions by mapping the full 
Anderson disorder model, 

H = -tY^ (ctcj + h.c) + Vim , (1) 

<i,j> i 

to an effective model of a cluster of Nc sites embedded in 
an effective (dynamical) bath with a complex self-energy 
(matrix of size N^). Here, we assume a binary alloy 
distribution for disorder, P{vi) = (1 — x)6{vi — va) + 
x5{vi — Vb), and further, that < ViVj > — < Vi >< 
Vj >= fij = C, a constant parameter. Strictly speaking, 
the SRO encoded in fij is a function of x, temperature 
and other variables depending on the specific physical 
situation under consideration, and in real materials, this 
dependence should be explicitly taken into account. 

In contrast to CPA, the method described below is 
tailor made to capture the dynamical effects of repeated 
scattering from a cluster of sites, which arc correlated in 
a manner described by fij over the cluster length scale. 
To proceed, we start with the embedded cluster Hamil- 
tonian, 

H = —t ^^(cgCc -|- h.c) + V Xana 

a a 

+ f^e^'=-«-(4cfc-|-M, (2) 

where 0, a = 1, ..,z denote a central site coupled (via 
t) to z nearest neighbors on a. d — z/2 dimensional lat- 
tice, and the last term describes the hybridization of the 
boundary of the chosen cluster with an effective medium 



(conduction electron bath fimction) that has to be self- 
consistently determined by a suitable imbedding proce- 
dure. We describe the details in course of the derivation 
below. 

Defining the diagonal and off-diagonal propagators on 
the cluster as Goo('*^) =< cqicJ >,Gao{'^) =< Ca',Co > 
,Gaa{^^) =< Ca',Ca >7 Start with the equation of 
motion (EOM) for G'oo(w): 

u}Goo{u}) = l + ty^^ Gao{u}) + v < xqCq; 4 > • (3) 

a 

Notice the appearance of a higher-order GF on the rhs 
of Eq. (3). It is naturally interpreted as the probability 
amplitude for having an electron at a site with disorder 
potential v. Its EOM reads, 

(iV-v) < XoCo; 4 > = < Xo> +t'Y< XoCa, 4 > • (4) 

a 

The EOM for Gao(w) on the rhs of Eq. (3) reads, 

wGao(w) = V< XaCa, 4 > +tGoo(w) + * X] ^J'O^'^) ' 

(5) 

and in a way similar to that leading to Eq. (4), we obtain, 

{cO — v) < XaCa\ 4 >= * < XaCQ\ c\ > +t ^ < XctCj;c\ > . 

(6) 

Continuing along with identical lines for the various 
higer-order GFs generated in Eqs. (3-5) give us 

U) < XoCa, 4 > = V < XoXaCa] cl > +t < XqCq; Cq > 

+ <xoCj;cl> , (7) 

U) < XaCo; 4 > = < Xa> +V < XaXQCo; 4 > +* < XaCa, 4 > ) 

(8) 

{ui -V) < XoXaCa, 4 > = * < XoXaCo] 4 > 

-I- i ^ < XoXaCj-, 4 > (9) 

and, 

{cO — v) < XaXoCo; 4 > = < XoXa > +t < XoXaCa, 4 > • 

(10) 

Finally, 

(W - Efe) < AoaCk; cl >= tk < AoaCa, 4 >; (H) 
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where Aoa = I, xq, Xa, xox^ for the various types of 
Green functions which couple the bath back to the cluster 
(see above), and < xoa >=< xqXc > is the non-local cor- 
relation function of the disorder potential over the cluster 
length scale (scaling like 1/d in d dimensions). For a sin- 
gle site cluster, we recover the exact CPA result using the 
EOM for Gao,Gko,< xoCo;cl > and < XoCfcicJ > only. 
Indeed, the local Green's function at the site is easily 
seen to be. 



1- < a;o > , <xq> 
Goo(w) = r-r-T- + 



ui — A{ui) uj — V — A(a;) ' 



(12) 



with A(a;) = Y.^, which is exactly the CPA 

result. 

In our cluster generalization, after a long and some- 
what tedious algebra, we finally find. 



Goo(w) = 



and. 



Gao(w) = 



l + V < XqCq;c\ > +{v/t)F2{uj) < XaCa;cl > 



UJ — zF2{(j) 



1 VOJ < XaCa, cj > +t(l + V < XqCq; cJ > 



(13) 



w - A(w) 



UJ — zF2{ij) 



(14) 



where 

V < xoco; >= 



v<xo>+{v/tr<xox^>^^^^ 



UJ — V — F2{uj) 



(15) 



and. 



with 



< XaCoc;c\ >= t ^i^i(a;) < XaCoicl > , (16) 



+ < Xa > — < XoXa > < XQXa > 

uj - Fi{uj) uj-v-Fi{uj) 



Here, Fiiuj) 



and F2(w) 



(17) 



-. Fi- 



uj-v-A{lj) """^ ^^y^! — uj-A(uj)- 

nally, the bath function, A{lo) is computed from the 
equation [11,12], 



r+w 



poie)de 



f 

^'^^^ l-w Goo^u;) + A(a;) - e 



(18) 



with A(u;) = t^Goo(w) for the Bethe lattice. At this 
point, one can show that the cluster-CPA technique de- 
veloped above always yields positive-definite local DOS, 
defined by p{uj) = —I'mGQQ{uj)/'K. To show this explic- 
itly, we observe that Goo(w) can be brought to a conve- 
nient mathematical form by simple algebraic manipula- 
tions 



Goo(w) = — — 

UJ - F2{UJ) 



<xa> 



+ 



UJ — V — F2{U)) 
< Xa > — < XoXa > , < XQXa > 

UJ - F2 (uj) 
< Xa > — < XfjXr. > 



UJ 



-Fi{u) 



■ + 



UJ ~ V ~ F2{uj) 

< XQXa > 

UJ — V — Fx{uj) 



(19) 



The first step in the derivation is to notice that each 
of the numerators is always positive definite by defini- 
tion. Clearly, to show that p{uj) is always positive defi- 
nite, we have now only to show that ImA{uj) < 0. From 

the EOM technique used above, A(w) = Y,k ZT^- ^ 



straightforward calculation shows that ImA{uj) < for 
any choice of the unperturbed DOS, po(e) > 0. Substi- 
tution in Goo(w) above immediately shows that the dis- 
order averaged DOS is always positive definite (clearly, 
self-consistency does not modify this conclusion). 

Few clarifications concerning the physical meaning of 
the set of equations is in order at this point. First, we 
notice that the carrier dynamics is an explicit function 
of the higher-order (in 1/d) SRO correlator, < xoXa >■ 
It is also easy to check that the system of equations are 
exact both in the band and the atomic limit, and the 
CPA result is readily recovered for the single-site cluster. 
To interpret the meaning of the bath function, A(w) in 
our approach, we begin by observing that one can view 
any selfconsistent cluster approximation as being a valid 
description in a regime with short-ranged order on the 
cluster length scale (analogous to the single site approx- 
imation being formally exact at mean field level). An 
exact solution of the problem implies consideration of an 
infinite cluster, and, of course, is an insoluble problem. 
Our choice for A(u;) above is then linked to the math- 
ematical consideration of short-ranged correlations over 
the cluster length scale only, or, equivalently, to the con- 
sideration of dynamical eff'ects of repeated scattering by 
a cluster consisting of a central site; plus its z nearest 
neighbors only. The effects of non-local SRO appears 
explicitly in the bath function A{uj) (i.c, in Eqn. (18)) 
via Goo{i^) as defined in Eqn. (13) with its explicit de- 
pendence on fij. It follows that the approach describes 
carrier dynamics in a situation where the carrier mean 
free path is of the order of the size of the chosen clus- 
ter (Z ~ a, the lattice constant for our cluster) in a fully 
selfconsistent way (see below), one step beyond the CPA 
where 1 = 0. 

It is interesting to notice that Goo{uj),Gao{i^) can be 
(formally) analytically expressed in terms of the corre- 
sponding diagonal and off-diagonal cluster sclf-cncrgics 
T,oo{uj),'£ao{i^) for d = 1,2,..., 00, as well as on certain 
special lattices. The above set of equations then consti- 
tute a closed set of simultaneous non-linear equations for 
the two self-energies, and are solved self-consistently to 
yield the renormalized (by disorder) DOS at the central 
site, p{uj) = — ^lmGoo(w). 
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The alloy correlation function (describing SRO) is 
given by < xoXa >=< xq >< Xa > +Coa in the gen- 
eral case, with Cqq encoding complete information about 
order-disorder instabilities in the alloy. It is important to 
notice that the dynamical effect of strong scattering by 
these short-ranged correlations (Coa) on the electronic 
self-energy is explicitly included within our formulation 
above. In particular, the electron can undergo repeated 
scattering on the atomic sites within the chosen cluster, 
and, depending upon the degree and character of SRO 
(see below), can be localized due to interference effects 
coming from repeated scattering from spatially separated 
centers; i.e, via Anderson localization. To address the is- 
sue of Anderson localization in our NLCPA scheme, we 
follow Economou et al. [10] and use the localization func- 
tion defined by: 



L{uj)^Kt I Goo(w) 



Goa{^)GaQ{^^) 



(20) 



where electronic eigenstates with energy u) satisfying: 

• L{lo) > 1 define extended states, 

• L{uo) < 1 define Anderson localized states, and, 

• L{uj) — 1 defines the mobility edge. 

Here, K is the connectivity of the lattice. The formalism 
developed above thus allows a complete determination 
of the Anderson transition and its dependence on lattice 
structure, type of SRO, and band-filling. 

III. NUMERICAL RESULTS AND DISCUSSION 



fact, the band split regime occurs for v/W > 1/4. How- 
ever, states near and at the band center become Ander- 
son localized before the band splits (Fig. 1) and the metal 
(incoherent)- Anderson insulator transition is continuous. 
For V < Vc, the incoherent metal has a very similar char- 
acter (breakdown of the quasiparticle) to that found in 
d — oo (CPA). It is also clear that the configuration aver- 
aged single particle DOS shows no anomalies across the 
Anderson localization transition, in agreement with well 
known [13] arguments. 




FIG. 1. Cluster-CPA DOS the binary-alloy distribution in 
the extreme SRO limit on a Bethe lattice for < xo >= 0.5, 
< xoa >~ 0.0 and various values of the local disorder po- 
tential, V = 0.1,0.2,0.3,0.4. Shaded regions define extended 
states, and unshaded regions define Anderson localized states. 



In this section, we describe the results obtained from 
the numerical solution of the self-consistent set of cou- 
pled nonlinear equations derived in the previous section. 
Since we are interested in generic effects of atomic SRO 
on carrier dynamics, we choose the semicircular unper- 
turbed DOS po{E) = :j^^VW'^ — E'^ as an approxima- 
tion to the actual DOS for a three dimensional cubic 
lattice [11]. This leads to a considerable simplification in 
the numerics without affecting the generic features qual- 
itatively. We work with W = 1.424 eV and study the 
fully renormalized DOS, p{ijj) and L{uj) as functions of 
the alloy composition y = (< xq > /(I— < xq >)), the 
atomic SRO parameter Cqo and the disorder strength v 
for a half-filled band. 

We begin with the symmetric case with y = \, and 
extreme random disorder, i.e, Cqq = or < xqq >=< 
Xq >< Xa >■ In d = cxD, this corresponds to the CPA, 
with the metal-insulator transition occuring continuously 
at V > W . Inclusion of SRO drastically changes the pic- 
ture. The M-I transition now occurs much earlier. In 
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FIG. 2. (Cluster-CPA DOS on a Bethe lattice for 

< XQa >— 0.0, V = 0.5 and different alloy concentrations, 

< xo >= 0.4,0.3, 0.2,0.1. 
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In Fig. 2, we show the effect of changing the alloy com- 
position on the critical value of v needed to split the band. 
For the case y = 2/3 (i.e, < xq >— 0.4), a larger Vc — 0.5 
is required to localize sates at Ep, and it increases to 
Vc = 0.85 for y = 3/7(x — 0.3). This Anderson insu- 
lating state (notice that the Fermi level, denoted by the 
vertical lines in our plots, lies in the region of localized 
states) is explicitly related to our inclusion of the effect of 
carrier scattering on short-ranged intersite atomic corre- 
lations Coq = C) and is never observed in the CPA solu- 
tion (d = oo), which always predicts an incoherent metal 
for a (half-filled band) particle-hole asymmetric disorder 
distribution. A continuous transition from an Anderson 
localized insulator to an incoherent metal is clearly seen 
upon decreasing y for a fixed disorder strength. 
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FIG. 3. Effect of partial SRO on the Cluster-CPA DOS 
for a Bethe lattice with < xq >= 0.5, v — 0.5 and 
< xoc >= 0.1,0.2,0.3,0.4. 

Next, we focus on the effect of varying the SRO pa- 
rameter on the electronic structure. In Fig. 3, we show 
the DOS for v = 0.5 for different values of < xoa >— 
0.1,0.2,0.3,0.4. Clearly, introducing partial order (ac- 
tually, this corresponds to increasing Coa) results in in- 
creased tendency to itinerance, moving Ep out of the 
region of localized states. The AL insulator to incoher- 
ent metal transition is clearly observed with increasing 
Cqq, and is a concrete illustration of an insulator- metal 
transition driven by the degree of atomic SRO (partial 
order) in a system. Clearly, increasing partial order (no- 
tice that increasing Coa corresponds to increasing the 
probability of having the same potential on the cluster 
sites) reduces the localizing effect of strong (repeated) in- 
tracluster disorder scattering, driving the I-M transition 
via increased itinerance. The situation is very analogous 
to the case where the pure Anderson disorder model is 
supplemented by additional short-range correlations in 



the hopping [14], where increasing the off-diagonal ran- 
domness drives an insulator-metal transition for a fixed 
diagonal disorder strength. 

Qualitatively similar behavior is seen for an asymmet- 
ric alloy distribution. In Fig. 4, we show the DOS for 
V — 0.5 and y = 3/7 (i.e, < xo >= 0.3). Interest- 
ingly, the spectrum shows additional features, but the 
AL insulator-metal transition with increasing Cqq follows 
the trend for the symmetric [y = 1) case. 
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FIC. 4. Effect of partial SRO on the Cluster-CPA DOS 
for a Bethe lattice with < xo >— 0.3, v — 0.5 and 
< xoa >= 0.0,0.1,0.2,0.3. 

Increasing the ratio v/W reveals rich structures in the 
DOS. In Fig. 5, we show the one-electron DOS for v = 1.5 
with < XQa >= 0.0,0.1,0.2,0.3. Next, we turn our at- 
tention to Fig. 6, which shows the evolution of the DOS 
for the asymmetric alloy distribution with < xq >— 0.3 
for the same parameters. In this case, we are already in 
the split-band regime. Very rich structure is seen in the 
results. For comparison, we know that the corresponding 
DOS obtained for these cases within CPA {d = oo, not 
shown) shows a split band structure with only upper- 
and lower "Hubbard" bands. Obviously, CPA is inca- 
pable of resolving the fine structure in the DOS originat- 
ing from repeated scattering between spatially separated 
scattering centers. The rich structures seen in the cluster 
generalization correspond partially to these effects, and 
can be traced back to the spectrum of eigenstates of the 
isolated cluster. In fact, the multiple sub-bands can eas- 
ily be shown to be centered around eigen-energies of the 
isolated cluster for the case of the symmetric alloy distri- 
bution with foa = 0. However (Figs. 5-6), in the general 
case with foa 7^ 0, one sees eight or nine distinct sub- 
band structures. We interpret the additional structures 
as arising from atomic SRO (non-zero foa)', in particular 
from "shake-up" effects originating from strong resonant 
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scattering of carriers (from the atomic SRO) from clus- 
ter sites. For comparison, we remark that couphng the 
two-site cluster to the "bath" (rest of the lattice) via sec- 
ond order processes in the hopping (corresponding to the 
"Hubbard I" approximation for the cluster) is incapable 
of accessing SRO effects in a consistent way. In particu- 
lar, in addition to violating the Hubbard sum rules [15], 
it cannot yield "shake-up" features in the DOS, always 
yielding only six bands centered around the eigenvalues 
of the 2-site cluster, and broadened by an amount 0{t). 
This discussion shows the importance of treating the ef- 
fects of both itinerance (via A{u})) and the (incoherent) 
resonant scattering on the same footing, and reveals the 
weaknesses inherent in uncontrolled approximations. 
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FIG. 5. Cluster-CPA DOS for < xo >= 0.5, v = 1.5 and 
< xoa >= 0.0,0.1,0.2,0.3. 
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FIG. 6. Cluster-CPA DOS for < xo >= 0.3, v = 1.5 and 
< xoc, >= 0.0,0.1,0.2,0.3. 



Additional interesting features observed from the cal- 
culations deserve comments. We clearly observe that the 
localization function L{uj) shows non-analytic behavior 
near each subband edge, but no non-analyticities across 
the mobility edge. More detailed characterization of the 
Anderson insulator-metal transition requires a detailed 
study of the two-particle responses [16] (density correla- 
tions, optical conductivity) and is left for future work. 

Our results show some resemblance to those obtained 
by Rowlands et al. [17] for the same model (in Id) using 
the Korringa-Kohn-Rostoker Non-Local Coherent Poten- 
tial Approximation (NLCPA) employing the dynamical 
cluster approximation (DC A). However, no attempt has 
been made there to study Anderson localization. Fur- 
ther, our treatment of partial SRO is very different from 
theirs. It is worth pointing out that our results, along 
with the NLCPA ones, are quite different from those ob- 
tained by Jarrell et al. [5]. However, we are presently 
unable to quantify the reasons behind these differences. 

We emphasize that the approach developed here has a 
wide applicability to various problems where the effect of 
atomic (chemical), magnetic, Jahn- Teller, etc., SRO on 
the character of carrier dynamics is an important issue. 
In particular, it should be applicable to the problem of 
electronic structure of disordered TM alloys [18], and to 
more recent cases such as hole-doped manganites [19], 
where a plethora of experimental work clearly demon- 
strates the importance of such effects in a correlated en- 
vironment. Applications to such systems is in progress 
and will be reported elsewhere. 
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